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We study the A I = 1/2 rule for kaon decays and the Bk parameter for K° — K° mixing in a 
dual 5-dimensional holographic QCD model. We perform, in the chiral limit, computations of the 
relevant four-point current-current correlators, which depend upon self-interactions among the 5D 
bulk fields. Spontaneous chiral symmetry breaking (ySB) i s realized through boundary conditions 
on the bulk fields. Numerical results are analyzed in comparison with QCD, chiral perturbation 
theory (x'PT) and data, finding reasonable agreement with the experimental values of the gg and 
<727 parameters describing the A7 = 1/2, 3/2 decay channels. 


I. INTRODUCTION 

The gauge/gravity duality 111 provides new tools 
to investigate non-perturbative QCD. Holographic dual 
models which allow the study of confinement and X'SB 
have been developed in essentially two different ap¬ 
proaches: models based on string theory/supergravity in 
10D Q, or phenomenologically inspired holographic dual 
models in 5D 0, §, Os, g. In the past year, a study 

of meson masses and decay constants in the 5D models 
have shown quite good agreement with data, at least for 
the leading meson excitations. These encouraging results 
are based on two-point current-current correlators, and 
in the 5D description the leading contributions to such 
correlators are obtained from bare 5D propagators. One 
can also compute n-point correlators at tree level; how¬ 
ever, the fact that 5D interactions must be taken into 
account makes an important difference that we shall in¬ 
vestigate here. The holographic calculation of four-point 
correlators is a non-trivial test for 5D dual models of non- 
perturbative QCD, with, perhaps, the most interesting 
aspect being the holographic computation of observables 
which are important in QCD, particularly the celebrated 
A I = 1/2 rule for kaon decays. Since our calculation is 
the first within this framework, we focus upon the sim¬ 
plest AdS/QCD model. 


II. THE A7 = 1/2 RULE, AND THE B K MIXING 
PARAMETER 

We start with a review of the relevant facts. Ne¬ 
glecting CP-violation there are two independent K° de¬ 
cays: 7\° —> 7T + 7r - and K° —> 7r°7r°. These two decays 
are combinations of A7 = 1/2 and A7 = 3/2 isospin 
amplitudes (Ao and Ai, respectively). Experimentally 
ReA 0 /ReA 2 = 22.2, and the unexpected largeness of 
this ratio is the A7 = 1/2 rule. In the chiral limit, 
and at 0 (p 2 ) in the chiral expansion, these two ampli¬ 
tudes are expressed in terms of the gs and g 2 7 parame¬ 
ters (see e.g. [§, 10)). In the following we will also dis¬ 
cuss the related but simpler observable, Bk, parameter¬ 
izing K° — K° mixing. Performing an OPE, its calcu¬ 


lation reduces to that of the hadronic matrix element 
< K°\Qas= 2 \K° > ( 7 *) = | fK M K B K {F) of the four- 
quark left handed current operator Qas =2 = 4?y y,uL^ d L^ d 

with L^ iq2 = gi 7 M ^ 1 ~ 2 75 ^ q -2 ■ Similarly, g s ,27 involve the 
K —> 7T7T matrix elements of two AS = 1 operators four- 
quark left handed current operators (above the charm 
threshold), Qi = 4 V^L^ d L^ u and Q 2 = Ar/^L^L^. 

The factorized contribution for these matrix elements 
is well known; the unfactorized contribution is more con¬ 
troversial. In the chiral limit the leading N c unfactorized 
contributions reduce to the calculation of integrals in Q 2 
(with Q the Euclidean momentum flowing between the 
left handed currents) of two distinct four-point current 
correlators EE]- For B k = B k (ji) ■ CasMf), 


b k(p) = \-J dQ 2 W lrlr (Q 2 )( M ), (1) 
W LRLR (Q 2 ) = J ^^(9). 


W LRL R (q) = jim * 3 J d A xd A yd A z jor-Wiv-*) 
(0\T{L^ d (x)Rj s (y)L^ d (0)Rj s (z)}\0)\ conn , (2) 

while gs and g 2 r are given by 


98(f) = zi(n) ( ABl l ^ - 1 ) + 


(3) 


8 Bk(f) 


^(1-^^- IdQ 2 "^7 ' M). 


> Wllrr(Q 2 ) 


15 


921 ( 9 ) = Ziin) +z 2 (y) 


4B k (f) 


(4) 


In the above, Cas= 2 (p) and 21 , 2 ( 7 *) are the Wilson coef¬ 
ficients of Qas ~2 and Q 1>2 , /?// q2 = qiY ^ 1 + 2 5 \ 2 , and 

W llrr (Q 2 ) J 

W llrr(9 ) = | im 0 * 3 j d A xd A yd A z e iqx+ik{ - v *) 
(0\T{L^ u (x)L I £ id (0)Rj u (y)Rd s (z)}\0)\ cona . (5) 
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III. THE HOLOGRAPHIC QCD MODEL 

To compute the four-point correlators through the 
gauge/gravity correspondence, we consider the Lq lQ2 and 
R ^ iq2 4D quark currents. By virtue of the AdS/CFT 
correspondence, these operators couple to the boundary 
values of 5D gauge fields L^(x v ,z) and R ,i (x l ',z), and 
the 5D gauge fields are massless. The 5D action is a 
SU(3)l x SU(3)r Yang-Mills theory 

~ j d^ d 4 x yfg {L a MN L a ’ MN + R a MN R°” MN ), (6) 

where M 5 is an undetermined mass scale, and M = (p,, 5). 
The fifth coordinate 2 runs from Lq to L\ 1 which are the 
positions of the UV and IR branes respectively. It is 
possible to set L 0 = 0 as in Refs. [3, IB] as our results 
are smooth as lo 0- We use the AdSs metric ds 2 = 
a { z ) 2 {r]f_iu dx^dx v — dz 2 ), where a(z) 2 = ( L/z ) 2 and r)p V 
is given by (+, —, —, —). The field strengths are Lmn = 
9mL n - d N L M - i[L M ,L N ], where L M = L a M T a , and 
similarly for Rmn- We normalise tr (T a T b ) = S ab , and 
define the vector and axial gauge bosons as Vm = ( Lm + 
Rm)/V2 and A M = ( L M — Rm)/V%, respectively. 

The holographic dual of this 5D theory is a 4D theory 
with three quark flavours and a global SU(3)l x SU(3)r 
chiral symmetry ill We work in the chiral limit 
and include the effect of the spontaneous ySB via the 
IR boundary condition on the axial vector fields II]. 
This is clearly a rather severe limit, yet the calculational 
simplification it affords makes it worthwhile. (The ef¬ 
fects of spontaneous and explicit ySB could be treated 
more precisely by the inclusion of a bulk scalar in the 
bi-fundamental representation of the gauge group SIo 


A. Propagators and 5D Interactions 

By the holographic correspondence, the UV bound¬ 
ary values of Vp and A M act as classical sources coupled 
to the 4D vector and axial global symmetry currents. 
Moreover, the 5D action evaluated on the solutions of 
the equations of motion (EOM) of the bulk fields defines 
the generating functional for current-current correlators 
in the 4D theory. We therefore need to solve the EOM 
with given UV boundary values of the fields, and sub¬ 
stitute back into the action. The EOM are non-linear, 
so we must solve them iteratively. This is the origin of 
Witten diagram construction for calculating the n-point 
Greens functions of operators in the dual theory I • The 
iterative solution of the EOM requires the bulk-to-bulk 
and bulk-to-boundary propagators for the gauge fields in 
the 5D theory. 

We work in the unitary gauge in both the SU ( 3)v and 
SU(3)a sectors. V 5 is eliminated completely in this gauge 
if we impose null Dirichlet boundary conditions for it 
on both branes. We impose warped Neumann boundary 
conditions on the A 5 field (i.e. d z {aA§) = 0 on both 


branes). In this gauge A 5 is a physical field because it has 
a zero mode which cannot be gauged away. The Vp and 
Ap propagators have both transverse and longitudinal 
parts in this gauge. Following Ref. E3 we have 

(WV") = —iGp ( 2 , z')P£ v - iG% ( 2 , z')P£\ (7) 

where P^A = rA v — and P£ v = The Ap 

propagator has a similar form, and both solve 




zSjz-z Q 
M 5 L 


( 8 ) 


The boundary conditions on ( 2 , z') are Dirichlet on 
the UV brane, Gp ( 2 , 2 ') | z=L = 0, and Neumann on the 
IR brane, d z G^{z,z')\, =L = 0. We impose Dirichlet 
boundary conditions on both branes in the axial sec¬ 
tor Gp(z,z') |_ Lo Li = 0, to encode ySB. The bulk-to- 

boundary propagator is defined as: 




(*') = -^d z (v»v v ) 


dADS 


2=io 


(with a similar equation for Ap), where 


(9) 


{V*V V ) 


(z') 

dADS 


—iKp (z’)P£ v - iK%(z')Pr, (10) 


K’ A (z') = -—d z GV’ A ( z ,z') z 


z=Ln 


( 11 ) 


and similarly for K^' A (z'). The solutions are given by: 


G V p A {z,z') z<z , = 7T22' X 
[Aj x [pz) + By^pz^GMpz') + pyiivz')} 
2M 5 L(AD - BC) 


,( 12 ) 


2 > z' implies the exchange 2 <-► z'. For the vector sec¬ 
tor the coefficients are A = — Vi(pTo), G = — Vo(pii), 
B = J\ (pLo) and D = J 0 (pLi); for Ap we obtain 
A = -yi(pL 0 ), C = -Vi(pii) B = JlipLo) and 
D = In the zero momentum limit: 


GV(z,z') z< z' =- 


2 M 5 L 


(z 2 -L 2 0 ), 


(13) 


(2 > z' again swaps 2 with z'). For the axial sector 


G A {z,z') z<z , 


1 (z 2 -L 2 )(z> 2 -Lj) 

2 M 5 L Ll - L\ 


(14) 


and swapping 2 and z' for 2 > z'. The bulk-to-boundary 
propagators are then 


<’V) 


2 ' CMpz') + Dy^pz') 

L 0 AD - BC 


(15) 


while Kq {z') = 1, and K A (z') = ■ The physi¬ 

cal condition of 4D current conservation eliminates the 
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FIG. 1: The three basic structures of the 36 individual 
Witten diagrams: X-diagrams involving the 4-boson vertex, 
Y-diagrams involving two 3-boson vertices, and Y diagrams 
where As propagates. The exterior circle represents the (UV) 
boundary of the AdSs space, the solid lines are boundary-to- 
bulk or bulk-to-bulk (axial or vector) gauge propagators, and 
the dotted line is a bulk-to-bulk As propagator. 


longitudinal piece of the bulk-to-boundary propagators. 
The positions of the poles of Kp and Kp are the masses 
of the vector and axial-vector mesons, respectively. This 
is the usual association of the masses of the Kaluza-Klein 
modes in 5D with the masses of the bound states of the 
dual theory 0 Note that this tree-level treatment of 
the 5D model incorrectly predicts that the mass of the 
n’th excited resonance scales like n (for n 1) as op¬ 
posed to the s/n Regge behaviour. This failing is to be 
expected, as the simple 5D theory here presented is in¬ 
valid in the far UV fl4 l. 

We also note that the current-current correlator in 
the axial sector as p —> 0 is IUi(p 2 )| p= o = F% = 
2 M 5 L/(L\ — Lq) -J, {B, ?6 |. The A 5 propagator in this 
gauge is ( A 5 A 5 ) = — iGp(z,z'), where Gp(z,z') is the 
limit as £ —> oo of the solution of 



£—dz + 2 + P 

z z 


G 5 p (z,z') 


zS(z-z') 

m 5 l 


(16) 

with the boundary conditions d z (aG p (z, z'))\ _ Lq l ^ = 0, 


giving 


G~{z, z') — 


—zz 


M 5 L{Lf - Ll) p 2 ' 


(17) 


The interactions are the three-gauge boson and four- 
gauge boson vertices derived from the 5D action Eq. ©• 


E x + Sa 5 + Sy, with 

S x = 3 iM 5 L J ^{[Kf + Kf][Kf + Kf ] 

~ Kq K p K p ) (18) 

Sa b = -|(M 5 L) 2 | ^^G 5 p (z,z')A'(z,z') (19) 


where A'(z, z') is a function whose explicit form is known. 
As for the Y-diagrams, the integrations are more in¬ 
volved, but the sum can be written as E y = Y\ + Yi + 
Y 3 + Y 4 - 2Y 5 - 2Y 6 , where 


Yi = 3 i 



dz d z' 

z z' 


Kp 1 {z)G b J (z, z')K? (z')At (z)I^ (z') (20) 


with { abcde} = {VVVVV}, {VAVAA}, {AAAVV}, 
{AVAAA},{AAVVA},{AVVAV} for i = 1,...,6. Per¬ 
forming the Wick rotation to Euclidean space, we obtain: 


Wlrlr(Q 2 ) = = iQ) = ~2 Wllrr(Q 2 )- (21) 

6 F4 

Thus we find that the ratio of the two correlators is inde¬ 
pendent of the integration in 2 -space, and is solely due to 
their different SU(3)l x SU(3)r symmetry factors. The 
integrals can be done analytically with finite values for 
L i and Lq, but our observables have a smooth Lq —^► 0 
limit, and much simpler expressions arise if we set L 0 = 0 
giving: 


E (Q) = 3iM 5 L 
69 


60 1 2 0 


+ 


32 


Q*L*I* 


32 16 69 

Q & L\ ~ 5 Q 4 L\ ~ mf 
3 ( 8 

15Q 2 L2/2 + 5 Q 4 Lf/ 0 2 
64 1 

Q e L\l 0 + QLGoh. 


+ 


2 

WWl 


( 22 ) 


where io,i = Iq,i(QLi) are modified Bessel functions of 
zeroth and first order respectively. 


B. Four-point current-current correlators 


In order to calculate the four-point current-current cor¬ 
relators we have to solve tree level four-point Witten dia¬ 
grams which contain either one four-gauge field vertex or 
two three-gauge field vertices. Note that certain 5d Wit¬ 
ten diagrams do not contribute to the four-point func¬ 
tions m and © because they do not preserve the flavor 
structure of the 4d theory. 

We find that the group theory structure is a common 
factor in all diagrams, so denoting by E the sum of all 
diagrams with group factors amputated we have E = 


IV. RESULTS AND DISCUSSION 


The three input parameters of our model are M$L : 
Lq and L\. From the expressions in Section II.A one 
can easily compute the boundary-to-boundary propaga¬ 
tor in the vector sector, and match its large Euclidean 
momentum behaviour to that known from perturbative 
QCD deriving the relation^, 5j 


M 5 L 


N c 

12t r 2 ' 


(23) 


One is then left with the two dimensionful parameters, 
1/Lq giving the UV cutoff scale, and 1/Li setting the 
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IR scale. As mentioned before, one problem of the sim¬ 
ple 5D holographic model presented is that it does not 
correctly reproduce the spectrum of meson resonances 
in the regime above approximately 1500 MeV. In addi¬ 
tion, in computing the physically relevant observables 
gs and (727 one must in practice choose a scale p at 
which to evaluate the Wilson coefficients z\(p) and Z 2 (/it). 
One simple choice we make in the following is to take 
1/Lq = 1500 MeV, imposing this as an upper hard cut¬ 
off when we evaluate the Q 2 integrals in Eqs. 0 - 0 , and 
identifying this hard cut-off with the short distance renor¬ 
malisation scale p. Although this matching procedure is 
quite crude, Refs. 0 have shown that such a proce¬ 
dure captures the dominant contributions. In a later 
longer publication we investigate in some detail other 
more sophisticated matching procedures. With respect 
to the \PT and perturbative QCD calculations of El 
and 0 , 0 of the Br and A/ = 1/2 rule observables, 
our model, if valid, presents the interesting possibility 
of calculating the intermediate Q 2 contribution around 
0.1-2GeV 2 directly, rather than fitting it from an inter¬ 
polation between the yPT low Q 2 regime [l5| and of the 
OPE high Q 2 regime 0, jTj,]. 

One immediate check on our calculation follows from 
the fact that the most general form for the Greens func¬ 
tion Wlrlr(Q 2 ) in the large N c limit can be written 
in terms of the masses, Mi, and residues of the meson 
resonances 00 


y- ( | Pi . 7i \ 

\(Q 2 + M 2 ) + (Q 2 + M 2 ) 2 + (Q 2 + M 2 Y) ’ 

(24) 

and similarly for Wllrr with residues a[, /?■ and 7 '. We 
have explicitly checked that our results for Wlrlr(Q 2 ) 
and Wllrr(Q 2 ) agree with this general form. 

We now turn to the comparison of our results with 
the experimental data. Given Eq. and our choice of 
Lq = 1/1500 MeV -1 , the one remaining free parameter 
of our model is L\. The data to which this must be fitted 
are, most importantly, the pion decay constant, F„, the p 
vector-meson mass, m p , the a\ axial-vector meson mass 
m ai , and the gs and (727 (or Br) parameters. In the low- 
Q 2 regime chiral perturbation theory calculations 00 
give the behaviour of the correlation functions 

Wlrlr(Q~) = 6 — 24(2/i + 5/2 + Z 3 + ^)~p2 "• 

” 7T 

W LLR r(Q 2 ) = -- + (—+ ■■■ ( 25 ) 

with U the standard chiral-Lagrangian coefficients. In 
our fit procedure we use the \PT results, Eci. ll25ll . in 
the integrals over Q 2 below 1/Ai, and the results of the 
holographic calculation, Eqs. ED and 11221 1 in the regime 
1/L 1 to 1/L 0 . 

In terms of M$L, L\ and Lq the predictions of the 5D 


holographic model for F w and m p and m ai 


are @ 


2 M 5 L 


(26) 


and in a good approximation in the range of interest 


2.12 (Li - 0.282L o ) 
L\ {Li — L q ) 
3.38 (Li - 0.085L o ) 
Ti (Li ~ L 0 ) 


(27) 

(28) 


(We use a precise interpolation to the exact result for our 
fits.) A non-trivial success of the 5D holographic QCD 
model is that in the Lq = 0 limit (to which 1/Lo = 1500 
MeV is a good approximation) the ratio m ai /m p ~ 1.6 is 
very well reproduced independently of the values of M$L 
and Li, as already noted in 11 

Numerically performing the Q 2 integrations in Eqs. CD, 
0 , and 0, as a function of 1/Li and utilizing, as appro¬ 
priate, the leading order values of the Wilson coefficients 
Z\ and (taken from 0 with A^ g = 325 MeV, and 
p = 1500 MeV) we can fit the final parameter Li to the 
full set of observables. The A I = 1/2,3/2 data requires 
gslobs = 5.1 and g 27 \obs_ = 0.29, although the values to be 
explained are modified to gs\obs = 3.3 and g 27 \obs = 0.23 
when one takes into account the enhancement already 
provided by the calculated 0(p‘ i ) chiral corrections Q. 
Fitting to gs and (727 and the values of F n , m p and m ai 
leads to the result 


L~ x « 302 MeV (29) 


which implies (in all cases normalized to the data for ease 
of understanding) 


m p \th 

77l p \obs 

9a\th 

98 I obs 


; 1 . 00 , 

0.49, 


77lai | th 
171 a 1 jobs 
927 \ th ^ 
927|obs 


« 1.04, 
1.39 


F\ L 

F„\ 


| th 
obs 


0.80, 

(30) 


and a value of Br (1500 MeV) « 0.54 (Br w 0.76 using 
the Wilson coefficient value 6^5=2(1500 MeV) = 1.42, 
see e.g. Ref. 0 ). Considering the relative crudity of 
the model and the use of the large- N c expansion of QCD 
this is a reasonable fit to the data and we find this result 
very encouraging. We emphasise that Li and L 0 are the 
only free parameters in this fit to the five observables of 
equation Eq. m . It is noteworthy that in the expected 
domain of validity of the model (i.e. in the Q 2 interme¬ 
diate region), the 5D interactions induce a large increase 
of gs and a suppression for ( 727 . 

On the other hand the holographic model predicts the 
low-Q 2 behaviour of Wlrlr and Wllrr to be identical 
in form to xPT but with disagreeing numerical values 


Wllrr 


Wlrlr 

2 


1 | M 5 L 29 Q 2 
4 + 4 24Tf + 


(31) 
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It is possible that this failing is due to the simplistic 
truncation of the AdS space in the IR. Alternatively, the 
fault might be with our approximations, particularly the 
treatment in the 5D model of ySB. 

It would be interesting to include a bi-fundamental 
bulk scalar field of mass-squared —3 JL 2 associated with 
the qq operator, as in the models 0, Q , and a massless 
singlet scalar associated with the glueball states gener¬ 
ated by the gluon field operator G^ V G ^ V , to see whether 
their inclusion leads to an even better fit to the full set 
of observables. 
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